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Abstract 

A synthesis of tlie present knowledge on gamma-ray emission from tlie magnetospiiere of a rapidly 
rotating neutron star is presented, focusing on the electrodynamics of particle accelerators. The combined 
curvature, synchrotron, and inverse-Compton emission from ultra-relativistic positrons and electrons, which 
are created by two-photon and/or one-photon pair creation processes, or emitted from the neutron-star 
surface, provide us with essential information on the properties of the accelerator — electric potential drop 
, along the magnetic field lines. A new accelerator model, which is a mixture of traditional inner-gap and 

■ outer gap models, is also proposed, by solving the Poisson equation for the electrostatic potential together 
' with the Boltzmann equations for particles and gamma-rays in the two-dimensional configuration and two- 
, dimensional momentum spaces. 

o 

:^ 

O I When a star has exhausted its nuclear fuel, the interior thermal pressure can no longer support its own weight. 

■ If the star is not very massive, it collapses to the point at which electron degeneracy pressure halts gravitational 
^ I collapse, and a white dwarf is formed. However, beyond the critical mass [1], even the electron degeneracy 

pressure becomes insufficient to support the star. In the collapsing stellar core, most of the electrons are 
absorbed by protons via inverse (3 decay and eventually the neutron degeneracy pressure halts the gravitational 
' collapse [2, 3, 4]. For a progenitor star to collapse directly to a neutron star, its mass should be more than 10 
H I times greater than our Sun. A neutron star, which can rotate with a period as fast as P ^ 1 ms, is considered 
to be the only candidate for the source of remarkably stable radio pulses, of which fastest period is 1.55 ms 
(PSR J1939-I-2134 [5]). When it was determined that the period of the Crab pulsar was slowly increasing [6], 
the identification of pulsars as rapidly rotating neutron stars was essentially confirmed [7, 8] . 

Having established pulsars as rotating neutron stars, the next thing to discuss is the pulsar characteristics 
deduced from observations. Since all known isolated pulsars show gradual spin down, there must be some 
mechanisms causing the neutron star to lose their rotational energy. The rotational energy-loss rate is given by 

E = mh^ ~{2TTflP/P^, (1) 

where the neutron star is rotating with angular frequency il = 2tt/ P and moment of inertia /. For typical 
high-energy pulsars, P ~ 0.1 s, P ~ 10~^'^ss~^, / ~ lO^'^gcm^; thus, we obtain E — —4 x lO'^^ ergs s~^, 
which is enough to explain the observed 7-ray luminosities, lO"^^'^ ergs s~^ < < 10'^^ ergs s^^. For an object 
as small as ~ 10 km in radius to experience such a large torque, it must have a strong coupling to their 
surroundings, most likely through magnetic fields intrinsic to the neutron star. If the neutron star is losing its 
rotational energy via magnetic dipole radiation, the spin-down luminosity is given by 

Ldip = fcf^VVc', (2) 

*Postal address: TIARA, Department of Physics, National Tsing Hua University, 101, Sec. 2, Kuang Fu Rd., Hsinchu, Taiwan 
300 



1 





1 400 40 

- 300 30 






200 40 
150 30 








: - 200 20 

- - 100 10 




1 


100 20 





I 


Crab 


Vela 


50 6 


Geminga 




- 30 

- 20 5 
10 


1 


II 


40 

4 

30 
20 

2 

10 


1 


1 



B1 706-44 B1951+32 B1 055-52 



Figure 1: Light curves of six 7-ray pulsars above 100 MeV (dotted lines, right-hand scale) and above 5 GeV 
(dark histogram, left-hand scale). Each panel shows one full rotation of the neutron star. From Thompson 
(2004). 

where k = 2 sin^ cti/S is commonly used for a vacuum rotating magnetic dipole with inclination ai with respect 
to the rotation axis. Equating Ldip with ~E, we can infer the magnetic moment, /i, of the star [9]. For 
P = 0.1 s, P = 10^^'^ s s~^, / = lO'^^gcm^, we obtain ^ = 2.6 x 10"^" G cm'' and the surface magnetic field 
strength, P, = 5.2 x 10^^(r»/10km)^'' G, as the Newtonian value at the magnetic pole. In what follows, we 
consider possible emission models of isolated, rotating neutron stars, after briefly mentioning the high-energy 
observations. 

2 Gamma-ray observations and physical processes 

The Energetic Gamma Ray Experiment Telescope and the imaging Compton Telescope aboard the Compton 
Gamma Ray Observatory have detected pulsed signals from at least seven rotation-powered pulsars (Crab 
pulsar [10, 11]; PSR B1509-58 [12]; Vela pulsar [13, 11]; PSR B1706-44 [14]; PSR B1951+32 [15]; Geminga 
pulsar [16, 11]; PSR B1055-52 [17]). Four of them (Crab, Vela, Geminga, and B1951-h32) exhibit double- 
peaked light curves[18, 19, 20] (fig. 1). Since interpreting 7-rays should be less ambiguous compared with 
reprocessed, non-thermal X-rays, the 7-ray pulsations observed from these objects are particularly important 
as a direct signature of basic non-thermal processes in pulsar magnetospheres. 

Spectral energy distributions (SEDs) offers the key to an understanding of the radiation processes. Ref. [17] 
compiled a useful set of broad-band SED for the seven 7-ray pulsars. The most striking feature of these i/F^, plots 
is the flux peak above 0.1 GeV with a turnover at several GeV. Various models [21, 22, 23, 24, 25, 26] conclude 
that these photons are emitted by the electrons or positrons accelerated above 5 TeV via curvature process. Such 
ultra-relativistic particles could also cause inverse-Compton (IC) scatterings. In particular, for young pulsars, 
their strong thermal X-rays emitted from the cooling neutron-star surface [27, 28, 29] efficiently illuminate the 
gap to be upscattered into several TeV energies. If their magnetospheric infrared (0.1-0.01 eV) photon field is 
not too strong, a pulsed IC component could be unabsorbed to be detected with future ground-based or space 
telescopes. 

The curvature-emitted photons have the typical energy of a few GeV. Close to the star (typically within a 
few stellar radii), such 7-rays can be absorbed by the strong magnetic field (> 10^^ G) to materialize as a pair. 
On the other hand, outside of this strong-field region, pairs are produced only by the photon-photon collisions 
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Figure 2: Sideview of pulsar magnetosphere. Left: two representative accelerator models in the open zone; 
center, global electric current due to the EMF exerted on the spinning neutron star surface when CI ■ fi > 0; 
right: current derived in the inner-slot-gap model (green arrow, see § 6.4). 

(e.g., between the curvature-GeV photons and surface- or magnetospheric- keV photons). The replenished 
charges partially screen the original acceleration field, . If the created particles pile up at the boundaries of 
the potential gap, they will quench the gap eventually. Nevertheless, if the created particles continue to migrate 
outside of the gap as a part of the global flows of charged particles, a steady charge-deficient region will be 
maintained. This is the basic idea of a particle acceleration zone in a pulsar magnetosphere. 

3 Representative emission models 

The pulsar magnetosphere can be divided into two zones: The closed zone filled with a dense plasma corotating 
with the star, and the open zone in which plasma flows along the open field lines to escape through the light 
cylinder (left panel of fig. 2). These two zones are separated by the the last-open magnetic field lines, which 
become parallel to the rotation axis at the light cylinder. Here, the light cylinder is defined as the surface where 
the corotational speed of a plasma would coincide with the speed of light, c, and hence its distance from the 
rotation axis is given by tulc = c/Q. If a plasma flows along the magnetic field line, causality requires that 
the plasma should migrate outward outside of the light cylinder. In all the pulsar emission models, particle 
acceleration takes place within the open zone. 

For an aligned rotator (i.e., sinai = 0), open zone occupies the magnetic colatitudes (measured from the 
magnetic axis) that is less than -^/r, / njLc ■ For an oblique rotator, even though the open-zone polar cap shape 
is distorted, 7r(r*/roLc) gives a good estimate of the polar cap area. On the spinning neutron star surface, from 
the magnetic pole to the rim of the polar cap, an electro-motive force, « ^"^B^fr^/c^ k, 10^^'^ V, is exerted. This 
strong EMF causes the magnetospheric currents that flow outwards in the lower latitudes and inwards along 
the magnetic axis (center panel in fig. 2). The return current is formed at large-distances where Poynting flux 
is converted into kinetic energy of particles or dissipated [30] . 

Attempts to model the particle accelerator have traditionally concentrated on two scenarios: Inner-gap 
models with emission altitudes of ~ lO^cm to several neutron star radii over a pulsar polar cap surface [31, 
32, 21, 33, 34] and outer-gap models with acceleration near the light cylinder [35, 36, 37, 38, 39]. It is worth 
noting that average location of energy loss should take place near the light cylinder so that the rotating neutron 
star may lose enough angular momentum [40], which indicates the existence of the inner gap must affect the 
electrodynamics in the outer magnetosphere. 

It is widely accepted from phenomenological studies that coherent radio pulsations are probably emit- 
ted from the inner gap that is located near the magnetic pole. Moreover, coherent radio pulsations pro- 
vide the primary channel for pulsar discovery (1560 pulsars as of April 2006, see ATNF pulsar catalogue, 
http://www.atnf.csiro.au/research/pulsar/psrcat/). However, it is, unfortunately, very difficult to reproduce 
them by considering plasma collective effects. Therefore, in this review, I will focus on incoherent high- energy, 
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non-thermal radiation, which is the product of magnetospheric gaps. 

It is worth mentioning other models than the inner and outer gap models. A charge depletion can alter- 
natively be a consequence of a stable charge void (i.e., no plasma region with B ■ E 0) 'm sl neutron-star 
magnetosphere [41, 42, 43]. Even though their works provides important hints on electrostatics, I do not intend 
to speculate on the dynamics of the formation of the void, to take account of magnetospheric currents, which 
are not considered in their works. There is another model that radiation takes place in a pulsar wind at roughly 
10 to 100 light cylinder radii from the star as a result of the magnetic energy dissipation [44, 45, 46], which is 
a reexamination of the non-axisymmetric striped pulsar wind [47, 48] . However, to discuss the reconnection in 
a relativistic plasma is beyond the scope of this brief review. 

4 Magnetospheric current determination 

We briefly discuss how to determine the magnetospheric current. Adopting the force- free limit (i.e., neglecting 
plasma inertia), Ref. [49] solved the Grad-Shafranov (GS) Eq. [50, 51, 52] {x^ - l){d^ + dl)P +{x+ l/x)d^P = 
jdj/dP, to obtain the magnetic flux function P, which determines the fieldline configuration, where x and 
z denote the distances from the rotation axis and the equatorial plane, respectively, in zulc unit, and j the 
magnetospheric current on the poloidal (i.e., meridional) plane. This work inaugurated numerous attempts [53, 
54, 55, 56] on solving the GS Eq. in recent several years, and k w 1+sin^ Oi is suggested in Eq. (2), by determining 
j{P)- However, it is noteworthy that they determined j by imposing a c;ontinuity of P at the turning point 
X = 1 (i.e., the light cylinder), and by simply neglecting the highest-oder term at a; = 1 (including a regular 
Taylor expansion, which gives only one of the two independent solution bases if we suppress the z dependence). 
Instead, we can adopt a singular perturbation theory to stretch the region a; « 1 by transforming P — Q/{x—l). 
Then there appears no singularity in the highest-order derivative term [57]. To find a particular solution, we 
can, for example, introduce a regulator such that P = Q/[x — 1 — e{z)], recovering the imaginary part by 
Fourier analysis from the equations of force-free electrodynamics [58, 55], where Q satisfies the GS Eq. without 
z derivatives and e is a small complex function. In this case, x = 1 + Re(e) will give the Alfvenic separatrix 
surface [59, 60, 61], after recovering particle inertia. Any way, we must connect the interior solution with the 
asymptotic ones in x < 1 and x > 1. It means that j, a free parameter in the force- free framework, cannot be 
constrained by the continuity of P at a; = 1, but should be determined by the global requirement including the 
dissipative region, which gives the electric load in the current circuit within the sub-fast-magnetosonic region 
(extended to large distances in a magnetically dominated magnetosphere) from causality requirement. 

5 Geometrical consideration 

In this section, we compare the emission geometry of a few representative gap models. For this purpose, we 
assume that the photons are outwardly emitted along the local magnetic field direction in the corotating frame 
of the magnetosphere. Then the lightcurve morphologies are subject to the general relativistic corrections (i.e., 
field distortion and light path bending near the star) and the special relativistic corrections (i.e., aberration of 
light, time-of-flight delays, and field retardation near the light cylinder). Aberration and time-of-flight delays 
leads to comparable phase shifts, ~ —w/wlc at distance w from the rotation axis. On the leading side (left 
panel of fig. 3), these phase shifts add up to spread photons emitted at various altitudes over 0.4 in phase. 
On the other hand, on the trailing side (right panel), photons emitted later at higher altitudes catch up with 
those emitted earlier at lower altitudes, arriving at an observer within a small phase range 0.1 and produce 
caustics [62, 39] in the phase plot and light curves [63, 64]. 

A single inner-gap beam, which is emitted from the lower altitudes, can produce a variety of pulse profiles. 
The top left panel in Fig. 4 represents the phase plot, — photon intensity in the direction specified by <I> (pulse- 
phase) and C (observer's viewing angle with respect to the rotation axis). An observer at C detects photon 
flux indicated by the shade at each pulse phase (bottom left panel). Because fades away near the perfectly 
conducting edge of the open zone, the gap is shorter near the pole and extends to higher altitude near the rim 
(right panel). Faint off-beam curvature radiation above the gap can be seen outside the main beam, as the faintly 
shaded regions indicate[65]. A great deal of effort has been made on the inner-gap model [66, 67, 68, 69, 70]; 
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Figure 3: Photon propagation direction in a rotating neutron-star magnetosphere (top view). Lejt. Leading 
side emission, which forms weak photon intensity; right: Traihng side emission, which forms a caustic — pihng 
up of photons emitted at different altitudes. 

however, one has to invoke a very small inclination angle to reproduce the wide separation of the two peaks 
(fig. 1). Thus, a high-altitude emission drew attention as an alternative possibility. 

Assuming that the gap extends from the stellar surface to the light cylinder and that the photon emissivity 
is uniform within the gap, Ref. [64] demonstrated the formation of the double peaks arising from a crossing two 
caustics, each of which is associated with a different magnetic pole (fig. 5). Subsequently, Ref. [71] examined 
polarization characteristics and found that fast swings of the position angle and minima of polarization degree 
can be qualitatively reproduced within their two-pole caustic model. This type of emission — slot gap emis- 
sion [72] — fills the whole sky and all phases in a lightcurve. Most observers will catch emission from the two 
poles if a\ > 30° (for 45° < ( < 125° in fig. 5). The dark features show the accumulation of photons because of 
the trailing side caustics (e.g., the thick curve at C < 130° and |<i>| < 50°), and because of the overlap between 
the trailing side of pole 1 and the leading side of pole 2 near the light cylinder (e.g., the thinner branch of the 
Y feature at 50° < C < 90° and 50° < $ < 110°). The main peaks come from the trailing side of each pole, 
interpeak emission from the leading sides. 

Geometrically speaking, an outer gap can be classified as a subset of the slot gap extended above the null 
surface. An observer can see, in this case, emission from the two outer gaps associated with only one pole. For 
example, in Fig. 6, if C = 80°, the dotted line shows that we observe the photons emitted from the two outer 
gaps (shaded regions in the right panel) associated with the pole 2. There is no emission outside the sharp 
peak edges; thus, the emission is invisible at C < 30° or C > 150° for any obliquity. The first peak originate 
near O.Otulc (upper shaded region in the right panel), while the second peak caustic is formed by the photons 
emitted from the trailing field lines (lower shaded region in the right panel). 

6 Electrodynamical consideration 

Next, let us consider the physical basis of the geometrical models discussed in § 5. To this aim, we have to 
derive the Poisson Eq. for the non-corotational potential, which is applicable to arbitrary gap models. 

6.1 Poisson Equation for Electrostatic Potential 

Around a rotating neutron star with mass M, the background space-time geometry is given by [73] 



ds^ = gttdt'^ + 2gt^dtdip + grrdr'^ + geedO'^ + g^^dip^ , 



(3) 



where 




(4) 




(5) 
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Figure 4: Phase plot, sample lightcurve, and a sketch of the accelerator location for the inner-gap model, for a 
typical inclination angle ai — 10°. the central zoom gives the gap extent relative to the star size, the dashed 
lines outline the null surface, the shading in the lightcurve and gap sketch is the same, the phase plot illustrates 
the change in lightcurve as seen by different observers and the aperture of the pulsed beams. From Grenier and 
Harding (2006). 




phase * 

Figure 5: Same figure as fig. 1 for the slot-gap model, for a typical inclination angle oi\ = 45°. From Grenier 
and Harding (2006). 
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Figure 6: Same figure as fig. 1 for the outer-gap model, for a typical inclination angle aj — 65°. From Grenier 
and Harding (2006). 



Tg = 2GM/c^ indicates the Schwarzschild radius, and a = in/{AIc) parameterizes the stellar angular momen- 
tum. At radial coordinate r, the inertial frame is dragged at angular frequency lu = —gtip/gipip = 0.15f2/45r6~^, 
where 145 = //lO"*^ erg cm^, and rg = r*/10km. 
Let us consider the Gauss's law, 



1 



Ait 
72"^' 



(6) 



where V denotes the covariant derivative, the Greek indices run over r, 0, ^J—g = \/ grrgeepi, = cr^ sin 9 and 
= g^^p — gttgtpip, p the real charge density. The electromagnetic fields observed by a distant static observer are 

given by [74, 75] Er = F^, Eg = Feu K = F^u = (gtt + fft^f^)Fe^/y=^, S« = {gtt+gtv,^)F^r/V^, = 

—p^Frel\f—g, where F^v = — 9^^^ and denotes the vector potential. 

Assuming that the electromagnetic fields are unchanged in the corotating frame, we can introduce the 

non-corotational potential \1/ such that 



,ip-nt), 



(7) 



where n = t,r, 9, (p. If FAt + ^Fa^p — {A — r,0) holds, is conserved along the field line. On the neutron- 
star surface, we impose Fgt + flFg^ — (perfect conductor) to find that the surface is equi-potential, that is, 
dg'i = dt'if + rJ^tp^* = holds. However, in a particle acceleration region, Fai + ^Fa^p deviates from and the 
magnetic field does not rigidly rotate. The deviation is expressed in terms of 4*, which gives the strength of the 
acceleration electric field that is measured by a distant static observer as 
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where the Latin index i runs over spatial coordinates r, 9, tp. 
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(8) 
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Substituting Eq. (7) into (6), we obtain the Poisson equation for the non-corotational potential, 



V g \ p-w / 



where the general relativistic Goldreich-Julian charge density is defined as 

^2 



Pgj 



47r\ 



-g 



pt. 



(9) 



(10) 



If p deviates from pcj in any region, E\\ is exerted along B. In the limit r fg, Eq. (10) reduces to the 
ordinary, special-relativistic expression [76, 77], 



Pgj 



f2 • B _^ (f2 X r) • (V X B) 



2ttc 



Aire 



(11) 



Instead of {r,6,(f), it is convenient to adopt the so-called 'magnetic coordinates" {s,6^,(p^,) such that s denotes 
the distance along the magnetic field line, 0* and the magnetic colatitude and the magnetic azimuth of the 
point where the field line intersects the stellar surface. Defining that ^* = corresponds to the magnetic axis 
and that = to the plane on which both the rotation and the magnetic axes reside, we obtain the following 
form of Poisson Eq., which can be applied to arbitrary magnetic field configurations [78]: 



1* 1 



'-g 



dx^ dx^ 



= g' 



dx^ dx^ 00 dx^ dx^ ko dx^ dx^ 



dr dr 



'-gg 

^2' 



dr 



+ g' 



n2 



86 do pIj dip dip 



'-gg 



dx^ 

'de 



i2^i' 



c g,pcp h) d^x 



pI pI 



dip"^ 



(12) 

(13) 

(14) 



where x = r, x = 0, x'^ = tp, x = ,s, x = 6^, and x"^ = ip^,. The light surface, a generalization of the light 
cylinder, is obtained by setting ko = gtt + '^gt^p^ + g^pip^^ to be zero [79, 80] Eq. (8) gives E\^\^ = —{d^/ds)e,,ip, ■ 
Magnetic field expansion effect [81, 82] is contained in g^*^' , g^''^' , gf-f- . 



6.2 Traditional Outer-gap Models 

Let us solve Eq. (12) for a vacimm case, p = 0, on (s,0») plane, neglecting Lp^ dependence. Supposing that 
copious charged particles freely migrate along the field lines outside of the gap, it is natural to consider that the 
lower and upper boundaries coincides with specific magnetic field lines. Considering the field lines that thread 
the stellar surface on = plane, we can express the lower and upper boundaries as 6** = ff^^^ and 0* = 
respectively. A magnetic field line is specified by the dimensionless quantity, h = {Q^'^'^ — 9>,)/9f^^] h = and 
h = hm = {9^'^^ — ^™™)/^™^ corresponds to the lower and upper boundaries. In the transversely thin limit 
(i.e., ft-m ^ 1), we solve the elliptic-type partial differential Eq. (12) by imposing = at h = 0, h = hm and 
s = 0, while ds'if = at an artificially chosen outer boundary (s = 1.4ti7LC, say), which does not affect the 
solution except for the vicinity of the light cylinder. We assume • /i > throughout this paper. 

As an example, we apply this method to the Crab pulsar, which has been studied from radio, optical. 
X-ray to 7-ray wavelength since its discovery [83, 84]. Its period and period derivative are P = 33.0 ms 
and P = 4.20 x lO'^^g g-i jf .^^g ^^^^^ = 2/8 in Eq. (2), the observed E = -4.46 x lO^^ergs s'^ gives 
/X = 3.80 x 10^°((i/2kpc)^. From soft X-ray observations, 180 eV is obtained[85] as the upper limit of the cooling 
neutron-star blackbody temperature, kT. 

Adopting kT = 100 eV, /z = 4.0 x 10^°G cm^, and the magnetic inclination ai = 70°, which is more or less 
close to the value (65°) suggested by a three-dimensional analysis in the traditional outer gap model (CHR) [63], 



8 



we obtain a nearly vacuum solution (fig. 7) for a geometrically thin case /im = 0.047. In the left panel, we present 
E^^{s,h) at five discrete heights. The solutions become similar to the vacuum one obtained in CHR. For one 
thing, the inner boundary is located slightly inside of the null surface. What is more, maximizes at the 
central height, h = and remains roughly constant in the entire region of the gap. The solved distributes 

almost symmetrically with respect to the central height; for example, the dashed and dash-dotted curves nearly 
overlap each other. The gap has no outer termination within the light cylinder. Since the inner boundary 
coincides with the place where 5* vanishes, the region between the star and the inner boundary has 4" > 0. 
Therefore, negative charges pulled from the stellar surface with = —dg'^ < populate only inside of the 
inner boundary, in which ^' < holds. Similar solutions are obtained for a thiner gap, < 0.047, even though 
there appears a small peak near the null surface, which is less important. 



distance= 2 kpc 

TTTT] 1 1 I I I lll| 1 1 I I I lll| 1 1 I I I lll| 1 1 I I I lll| 1 1 I I I I II 



incl. angle = 70degJ 




distance along field line photon energy [eV] 



Figure 7: Traditional outer-gap solution obtained for the Crab pulsar with ai — 70° and — 0.047. Left: 
The field-ahgned electric field at discrete heights h ranging from 2h^/16, 5/im/16, Shm/IG, Ilhm/IG, 14/ii„/16, 
with dashed, dotted, solid, dash-dot-dot-dot, and dash-dotted curves, respectively. The abscissa indicates the 
distance along the field line from the star in the unit of the light-cylinder radius. The null surface position at 
the height h — /im/2 is indicated by the down arrow. Right: Calculated phase-averaged spectra of the pulsed, 
outward-directed 7-rays. The flux is averaged over the meridional emission angles between 44° and 58° (solid) , 
58° and 72° (dashed), 72° and 86° (dash-dotted), 86° and 100° (dotted), 100° and 114° (dash-dot-dot-dotted), 
from the magnetic axis on the plane in which both the rotational and magnetic axes reside. 

Refs. [35, 36] first considered this kind of vacuum solutions and suggested two outer gaps can be formed 
for each magnetic pole (i.e., the upper and lower shaded regions in the right panel of fig. 6). Considering a 
fan beam (instead of a pencil beam in the inner-gap model or a funnel beam in the inner-slot-gap model) as 
the emission morphology, they discussed the formation of cusped photon peak as shown in Fig. 6. Extending 
this morphological emission model, Stanford group [39] discussed the observed properties of individual 7-ray 
pulsars, their radio to 7-ray pulse offsets, and the radio- vs. 7-ray detection probabilities. Assuming (x 
and a power-law energy distribution accelerated e^'^'s, Ref. [22] estimated the evolution of high-energy flux 
efficiencies and beaming fractions to discuss the detection statistics. Another group in Hong Kong [23] examined 
the minimum trans-field thickness, hm, of the gap, imposing that the 7-7 pair creation criterion is met. They 
estimated the soft photon field emitted from the heated polar-cap surface by the bombardment of gap-accelerated 
charged particles, adopting essentially the same £'|| solution as Fig. 7. Extending this work, Ref. [63] developed 
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a three-dimensional outer magnetospheric gap model to examine the double-peak light curves with strong inter- 
pulse emission (fig. 6), and estimated phase-resolved 7-ray spectra by assuming that the charged particles are 
accelerated to the Lorentz factors at which the curvature radiation-reaction force balances with the electrostatic 
acceleration. 

The outer gap models of these two groups have been successful in explaining the observed light curves, 
particularly in reproducing the wide separation of the two peaks (fig. 1), without invoking a very small inclination 
angle (as in inner-gap models). However, if we solve Eq. (12) self-consistcntly with the particle and 7-ray 
Boltzmann Eqs.[78], we find that the 7-ray flux obtained for hm < 0.047 (i.e., traditional outer-gap models) is 
insufficient (right panel of fig. 7). Thus, we have to consider a transversely thicker gap, which exerts a larger 
i^ll because of the less-efficient screening due to the two zero-potential walls at ft = and /im- If the created 
current increases due to the increased /im, the gap inner boundary deviates the null surface and shifts inwards 
to touch the stellar surface at last[86]. Therefore, we will analytically prove this behavior in § 6.3 then move on 
to the main task. 

6.3 Inner boundary position 

In the transversely thin (i.e., vacuum) limit, 9» derivatives dominate in Eq. (12) to give —{B/B^)r~^d'^^/d9'^ « 

—AnpGj- Noting that pcj > and dspcij > hold within the gap, and that = at both 6^, = 6*™'" and 
0* = O'^^^, we find \1/ < and i^y = —dg^ > 0. That is, the negativity of p — pcj results in a positive E^^ in 
the gap, if ^, derivatives dominate. Then, what happens near the inner boundary where s derivatives become 
important. To avoid the reversal of the sign of E^^, we find that dE\\/ds = —d^'^/ds^ w —Anpcj > must 
be satisfied. Therefore, we can conclude that the inner boundary should be located slightly inside of the null 
surface where pcj changes sign. 

The same argument holds for a non-vacuum gap. At the inner boundary, dgE^^ w 47r(p — pcj) > leads to 

PGJ ^ uj\B, \p\ O 

-^^2^cV-W^^^^2^c' ^^^^ 

so that may have a single sign. It follows that the polar cap, where ^ B* holds, is the only place for 
the inner boundary of the 'outer' gap to be located, if the created particle number density is comparable to the 
Goldreich- Julian value at the surface, namely, p flB*/{2Trc). Such a non- vacuum gap must extend from the 
polar cap surface (not from the null surface as traditionally assumed) to the outer magnctosphere. On these 
grounds, we have to merge the inner- and outer-gap models, which have been separately considered so far. 

6.4 Inner-slot-gap models 

In traditional inner-gap models, the predicted beam size of radiation emitted near the stellar surface is too 
small to produce the wide pulse profiles (fig. 1). Seeking the possibility of a wide hollow cone of high-energy 
radiation due to the flaring of field lines, Ref. [87] first examined the particle acceleration at the high altitudes 
along the last open field line. This type of accelerator, or the slot gap, forms because the pair formation front, 
which screens E^^ , occurs at increasingly higher altitude as the magnetic colatitude approaches the edge of the 
open field region [88]. MHOS [89] extended this argument by including two new features: acceleration due to 
space-time dragging, and the additional decrease of at the edge of the gap due to the narrowness of the slot 
gap. Moreover, to give the physical basis of (purely geometrical) two-pole caustic model (§ 5), MH04a,b [90, 91] 
matched the high-altitude slot gap solution to the solution obtained at lower altitudes by MHOS, and found 
that the residual E^^ is small and constant, but still large enough at all altitudes to maintain the energies of 
electrons, which are extracted from the star, above 5 TeV. 

We should notice here that their inner-slot gap model is an extension of the inner-gap model into the outer 
magnetosphere, assuming a spacc-cliargc-limitcd flow (SCLF) that the plasma flowing in the gap consists of 
only the charges extracted from the stellar surface. This assumption is self-consistently satisfied in their model, 
because pair creation in the extended slot gap occurs at a reduced rate. However, the SCLF leads to the electric 
current that is opposite to the global current flow patterns, if the gap exists near the last-open field line (right 
panel of fig. 2). 
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6.5 Modern outer-gap models — super Goldreich-Julian current with ion emission 
from stellar surface 

We are, therefore, motivated by the need to contrive an accelerator model that predicts a consistent current 
direction with the global requirement. To this aim, it is straightforward to extend recent outer-gap models, 
which predict opposite to polar-cap models, into the inner magnetosphere. Extending the one-dimensional 
studies [92, 93, 94, 95, 96, 86], Refs. [97, 98] solved (a simpler form of) Eq. (12) to reveal that the gap inner 
boundary is located inside of the null surface owing to the pair creation within the gap, assuming that the 
particle motion immediately saturates in the balance between electric and radiation-reaction forces. 

Extending Ref. [97, 98] by solving the Lorentz- factor and pitch-angle evolution of individual e^'s, H06 [78] 
solved the gap electrodynamics for the Crab pulsar. For /im > 0.047, the created current density je becomes 
super Goldreich-Julian in the sense that p < pcj < holds at the inner boundary (left panel of fig. 8). The 
predicted 7-ray flux is much larger than the sub-GJ case of hm < 0.047 (i.e., traditional outer-gap solution). 
Moreover, the copious pair creation leads to a substantial screening of E^^ in the inner region, as the right panel 
shows. In this screening region, p{s, h) distributes (fig. 9) so that may virtually vanish. Because of the ion 
emission from the stellar surface, the total charge density p is given by p = pe + Pion, where pe denotes the sum 
of positronic and electronic charge densities, while pion does the ionic one. We should notice here that even if 
E\\ « occurs by the discharge of created pairs in most portions of the gap, the negative po — pcj inevitably 
exerts a strong positive £^|| at the surface (see H06 for details), thereby extracting ions from the surface for the 
solution with super-GJ current. On these grounds, we can regard this modern outer-gap solution as a mixture 
of the traditional inner-gap model, which extracts electrons from the surface with i^y < 0, and the traditional 
outer-gap model, which exerts positive E\\ because of the negativity of p — pgj- 

Most of the inward-directed 7-rays materialize outside of the gap, leading to a pair cascade by colliding 
with the surface X-rays. Theses pairs quickly lose initial inward momenta by a strong synchrotron radiation to 
become non-relativistic, then are resonantly scattered by the surface X-rays, and finally form a pulsar wind. It 
follows that 8 X 10"^* particles are created per second, which appears less than the constraints that arise from 
considc^ration of magnetic dissipation in the wind zone[99] (10''° s~^), and of Crab Nebula's radio synchrotron 
emission[100] (lO^"""^ s^i). 

I summarize accelerator models/theories in Table 1. The second and third columns represent the spatial 
dimension of the model and the sign of exerted E'n . 
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Figure 8: Modem outer-gap solution obtained for the Crab pulsar with ai = 70° and > 0.048: The dotted, 
solid, dashed, dash-dotted, and dash-dot-dot-dotted curves corresponds to — 0.047, 0.048, 0.060, 0.100, and 
0.160, respectively. Left: Created current density jc (in unit of Q.B /2tt) as a function of the transfield thickness 
h. If jo appears below (or above) the dashed line, \cpG.-i / [VlB l2i:)\s=o, the solution is sub- (or super-) GJ 
current. Right: (s, /i,„/2). 

6.6 Inner-slot-gap model vs. modern outer-gap model 

It is worth comparing the results of II06 with the inner-slot-gap model proposed by MHOS, MH04a,b, who 
obtained a quite different solution (e.g., negative E\\ in the gap). The only difference is, essentially, the transfield 
thickness of the gap. Estimating the transfield thickness to be A^sg ~ hraf* \J r/wi^Ci which is a few hundred 
times thinner than II06, they extended the solution near the polar cap surface [82] into the higher altitudes 
(towards the light cylinder). Because of this very small A^sg, emitted 7-rays do not cfhciently materialize 
within the gap; as a result, the created and returned positrons from the higher altitudes do not break down the 
original assumption of SCLF near the stellar surface. 

To avoid the reversal of E\\ in the gap (from negative near the star to positive in the outer magnetosphere) , 
or equivalently, to avoid the reversal of the sign of the effective charge density, poff = p — PGii along the field line, 
MH04a,b assumed that Pcs/B nearly vanishes and remains constant above a few neutron star radii. Because 
of this assumption, is suppress at a very small value, which reduces pair creation significantly. To justify 
this p/ B distribution, they proposed that p should grow by the cross field motion of charges due to toroidal 
forces [101], and that Pcs/B is a small constant so that cpcs/B may not exceed the flux of the emitted charges 
from the star, which ensures the equipotentiality of the slot-gap boundaries. Because of the very thin gap, 
A^sG ~ the cross-field motion becomes, indeed, important. However, the assumption that Pcs/B is a 

small positive constant may be too strong, because it is only a sufficient condition of equipotentiality. 

On the other hand, H06 considered a thicker gap, A/sg ~ 0.5/i„i1I7lc, in which the cross-field motion 
is negligible. In the inner magnetosphere, Pcs/B becomes approximately a negative constant, owing to the 
discharge of the copiously created pairs; as a result, a positive is exerted. For a super-GJ solution, we obtain 
jc +iion ~ 0.9 > pcft/ {i^B /2'!r), which guarantees the equipotentiality. 

In short, whether the gap solution becomes MII04a way (with a negative as an outward extension of 
the inner-gap model) or II06 way (with a positive as an inward extension of the outer-gap model) entirely 
depends on the transfield thickness and on the resultant Pcs/B variation along the field lines. If AZsg ^ r^/10 
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Figure 9: Modern outer-gap solution of the total (solid), created (dash-dotted), and Goldreich- Julian (dashed) 
charge densities in ilB{s, h)/{2-Kc) unit, for ai = 70° and = 0.060 at three transfield heights, h. Because of 
an ion emission from the stellar surface, the total charge density deviates from the created one. 

holds in the outer magnetosphere, Pcs/ B could be a small positive constant by the cross- field motion of charges 
(without pair creation); in this case, the current is slightly sub-GJ with electron emission from the neutron 
star surface, as MH04a,b suggested. On the other hand, if A/sg > ci7Lc/40 holds in the outer magnetosphere, 
Pcff/B takes a small negative value by the discharge of the created pairs (fig. 9); in this case, the current is 
super-GJ with ion emission from the surface, as demonstrated by H06. Since no studies have ever successfully 
constrained the gap transfield thickness, there is room for further investigation on this issue. 
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